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ABSTRACT 

Context. The magnetic reconnection in the solar corona results in an impulsively gen¬ 
erated Alfven waves which drive a polar jet. 

Aims. Using the Hinode/EIS 2” spectroscopic observations, we study the intensity, ve¬ 
locity, and FWHM variations of the strongest Fe XII 195.12 A line along the jet to 
find the signature of Alfven waves. We simulate numerically the impulsively generated 
Alfven waves within the vertical Harris current-sheet, forming the jet plasma flows, 
and mimicking their observational signatures. 

Methods. Using the FLASH code and the atmospheric model with embedded weakly 
expanding magnetic field configuration within a vertical Harris current-sheet, we solve 
the two and half-dimensional (2.5-D) ideal magnetohydrodynamic (MHD) equations 
to study the evolution of Alfven waves and vertical flows forming the plasma jet. 
Results. At a height of ~ 5 Mm from the base of the jet, the red-shifted velocity com¬ 
ponent of Fe XH 195.12 A line attains its maximum (5 kms“^) which converts into 
a blue-shifted one between the altitude of 5 — 10 Mm. The spectral intensity conti- 
nously increases up to 10 Mm, while FWHM still exhibits the low values with almost 
constant trend. This indicates that the reconnection point within the jet’s magnetic 
field topology lies in the corona 5 — 10 Mm from its footpoint anchored in the Sun’s 
surface. Beyond this height, FWHM shows a growing trend. This may be the signature 
of Alfven waves that impulsively evolve due to reconnection and propagate along the 
jet. From our numerical data, we evaluate space- and time- averaged Alfven waves 
velocity amplitudes at different heights in the jet’s current-sheet, which contribute to 
the non-thermal motions and spectral line broadening. The synthetic width of Fe XH 
195.12 A line exhibits similar trend of increment as in the observational data, possibly 
proving the existence of impulsively generated (by reconnection) Alfven waves which 
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propagate along the jet. 

Conclusions. The numerical simulations show that the impulsive perturbations in the 
transversal component of velocity at the reconnection point can excite the Alfven waves. 

These waves can power the plasma jet higher into the polar coronal hole as vertical 
plasma flows are also associated with these waves due to pondermotive force. The sim¬ 
ulated Alfven waves match well with the observed non-thermal broadening along the 
jet, which may provide direct spectroscopic evidence of the impulsively excited Alfven 
waves within the polar jet. 

Key words. Magnetohydrodynamics (MHD) - Sun: atmosphere - Sun: corona - 
Methods: numerical 

1. Introduction 

Polar coronal jets are well observed, large-scale, and confined plasma transients in the 
solar atmosphere (see Nistico et al. 2009; Nistico & Zimbardo 2012). They can signifi¬ 
cantly contribute to the energy transport and the formation of nascent supersonic wind. 
The exact driving mechanisms of such jets are still debated. Broadly speaking, the two 
candidates, namely magnetic reconnection and magnetohydrodynamic (MHD) waves, are 
known to trigger these jets (e.g., Shibata 1982; Yokoyama & Shibata 1995; Cirtain et al. 
2007; Nishizuka et al. 2008; Filippov et al. 2009; Pariat et al. 2009, 2015, 2015; Srivastava 
& Murawski 2011; Kayshap et al. 2013, and references cited therein). The observations 
of Alfven waves associated with the polar coronal jets are either related to the imaging 
observations of photospherically driven these waves propagating along the jet (Cirtain et 
al. 2007) or some transversal perturbations evolved during the reconnection at the base 
of the jet (Nishizuka et al. 2008). Such Alfven waves are also ubiquitous in the 
localized magnetic structures (e.g., spicules, prominences, small-scale chromo¬ 
spheric fluxtubes) as well as in the large-scale corona (e.g., De Pontieu et al. 
2007; Okamoto et al. 2007; Tomczyk et al. 2007, Jess et al. 2009; Mathioudakis 
et al. 2013, and references cited therein). 

Detection of Alfven wave in solar magnetic structures is not yet well established as these 
waves are incompressible. In the polar coronal plasma, these waves are observed in the form 
of spectral-line profile variations (Banerjee et al. 1998, Harrison et al. 2002). Narrowing 
of the spectral line-width is attributed to the dissipation of small-amplitude Alfven waves 
(Harrison et al. 2002; O’Shea et al. 2005; Bempord & Abbo 2012; Dwivedi et al. 2014, and 
references cited therein). The growth and dissipation of Alfven waves have been modelled in 
the corona, which also explain the observed line-width variations (e.g., Pekenulu et al. 2002; 
Dwivedi & Srivastava 2006; Chmielewski et al. 2013, 2014). In particular, Chmielewski et 
al. (2013) have reported that the observed line broadening, particularly in the polar corona, 
can be explained in terms of impulsively generated non-linear Alfven waves. Non-linear 
Alfven waves may be the likely candidates for transporting energy in the solar corona 
(Murawski et al. 2015a). These waves can also power the large-scale, as well as confined 
plasma transients, e.g., solar jets, nascent wind, and plasma flows (Murawski et al., 
2015b). 
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Fig. 1. The intensity, Doppler velocity, and FWHM maps of Fe XII 195.12 A line, showing 
the polar jet moving off the limb outward. The white line is the position of the slit along 
the jet upon which the spectral parameters are estimated. The footpoint of the jet where 
slit’s lowest part lies, has the coordinates of (—28”, 964”). The white rectangular box is 
chosen at the limb of north polar corona from where the integrated line-profile of Fe XII 
195.12 A is derived and fitted in order to obtain its rest wavelength. 


It has been reported that the Alfven waves power various jets (e.g.. X-ray jets, spicules) 
and exist in a variety of coronal magnetic structures. However, they are identified as kink 
waves (Van Doorsselaere et al. 2008a, 2008b). The existence of mixed radial and azimuthal 
waves are also reported in theory and observations (Goossens et al. 2009, 2012; Tian et al. 
2012; Srivastava & Goossens 2013; and references cited therein). Numerical simulations 
support the general interpretation of the observed oscillations as a coupling of 
the kink and Alfven waves (see e.g. Pascoe et al. 2010, 2011). Kamio et al. (2010) 
have reported the spectroscopic observations of a rotating coronal jet and interpreted it 
as an evolution of kink waves and instability. The solar jets are found to be driven by 
various physical processes, e.g., the direct magnetic reconnection generated Lorentz 
force (Nishizuka et al. 2008), the reconnection generated pulse (Srivastava & Murawski 
2011), emergence and internal reconnection in small-scale kinked flux-tubes (Kayshap et 
al. 2013), etc. However, the pure Alfven waves driven coronal jets are difficult to detect 
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due to observational constraints and the dynamic nature of the jet’s typical magnetic field 
and plasma configuration. 

In the present paper, we study the Hinode/EIS spectroscopic observations 
of a polar jet. We find the spectroscopic signatures of impulsively generated 
Alfven waves and associated plasma flows within the jet. We model the ob¬ 
served physical processes (waves and flows) in the jet as a natural consequence 
of reconnection generated velocity pulse in the vertical and gravitationally 
stratified Harris current-sheet lying in the appropriate model atmosphere with 
VAL-III C temperature (Vernazza et al. 1981). The structure of the present paper 
is as follows. In Sect. 2, we present the observational results. Sect. 3 describes numerical 
model, governing equations, initial conditions, perturbations, and numerical solutions. In 
Sect. 4, the numerical results are shown. Discussion and conclusions are outlined in the 
last section. 

2. Hinode/EIS Observations of a Polar Coronal Jet 

2.1. Spectroscopic data and observational analyses 

A polar jet is observed using 2” slit scan by EUV Imaging Spectrometer (Culhane et 
al. 2006) onboard Hinode on 22 April, 2009. It is an established fact that the 40” and 
266” slots observe the temporal image data, while 1” and 2” slits are appropriate for 
the spectroscopic observations of the solar corona and Transition Region (TR). The EIS 
observes high-resolution spectra in two wavelength intervals, i.e., 170 — 211 and 246 — 
292 A using its short-wavelength (SW) and long-wavelength (LW) CCDs respectively. The 
observed data contain spectral line profiles of Fe XIII 202.04 A, Fe XI 188.23 A, Fe XV 
284.16 A, Fe XII 195.12 A, Fe X 184.54 A, Si VII 275.35 A, O V 192.9 A and He II 
256.32 A. The scanning started at 05:31:43 UT and ended at 06:33:29 UT. The exposure 
time on each scanning step was 31 s. The polar jet and related plasma column that was 
moving off the limb is scanned fully in a single spatio-temporal step at 06:29 UT because 
it went straight into the corona. This provides us an opportunity to understand the wave 
activity along the jet that already reached up to a certain height in the polar corona. The 
2” slit started scanning steps over the polar coronal hole with (Xcen, Ycen) ^ (30”, 939”). 
The observation window on the CCDs is 400 pixels high along the slit with a width of 40 
pixels in the horizontal direction. The Y-direction covers the solar atmosphere from 740” 
to 1139” (400 pixels with l”/pxiel). The direction of dispersion has the spectral resolution 
of 0.02 A. 

We apply standard EIS data-reduction procedures and calibration files to the data ac¬ 
quired at the EUV-telescope which is the raw (zeroth-level) data. The subroutines can be 
found in the sswidl software trecj^. These standard subroutines correct for dark-current sub¬ 
traction, cosmic-ray removal, flat-field correction, hot pixels, warm pixels and bad/missing 
pixels. The data are saved in the level-1 data file, while associated errors are saved in the 
error file. We choose the strong line Fe XII 195.12 A to examine the spatial variations of 
the intensity, Doppler velocity, and FWHM along the jet from its base in order to get the 
^ http://www.darts.isas.jaxa.jp/pub/solar/ssw/hinode/eis/ 
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clues of transversal waves. We co-align the Fe XII 195.12 A map w.r.t. the long-wavelength 
CCD observations of He II 256.32 A by considering it as a reference image and by estimat¬ 
ing the offset. The orbital and slit-tilt are also corrected for the data using the standard 
method described in the EIS software notes. We perform the double Gaussian fitting for 
the removal of the weak blend of Fe XII 195.18 A line that affects the line profile of Fe XII 
195.12 A. The fitting function is a Gaussian for line profile (see gauss_.pro in SolarSoft), 
while straight-line for background continuum (see line_.pro in SolarSoft). The fitting of 
the observed spectral line profile gives peak intensity, centroid (measure of 
flows), and Gaussian width (measure of thermal and non-thermal motions). 
We apply the procedure described by Young et al. (2009) in this context, which is also 
available in the EIS Software Note 17. We constrain Fe XII 195.18 Aline to have the same 
width as of Fe XII 195.12 A line. We search for the blend line Fe XII 195.18 A within the 
range of -1-0.06 A w.r.t. the centeroid of the main line Fe XII 195.12 A. We assume and 
fix that the contribution of the blended line will be searched within maximum 28 % limit 
compared to the peak intensity of the main Fe XII 195.12 A line. We perform the fitting 
over binned data of 2 pixel x 6 pixel to increase the signal-to-noise ratio and to obtain the 
resonable fitting and estimated parameters (cf., Fig. 1 and Fig. 2). The example of the Fe 
XII 195.12 A fitted spectra at 14 points along the chosen slit along jet are shown in Fig. 2. 
These fitted profiles of Fe XII 195.12 A are free from any contribution of the weak blend 
of Fe XII 195.18 A. The 14 spatial points along the jet corresspond to various heights for 
which the spectral parameters are derived (Fig. 3). 

It should be noted that the wavelength calibration and estimation of the reference 
wavelength of Fe XII 195.12 A line is performed using the limb method (Peter & Judge 
1999). The box is chosen at the north pole limb from where the integrated spectra is 
derived (cf., Fig.l, intensity map). It is likely that in and outward motions near the limb 
will almost cancel to each other. Therefore, the centroid of the line profile will represent 
the rest wavelength of particular emission line. The Gaussian fitting gives the estimate of 
the reference wavelength as 195.128 A. This reference wavelength is used in deriving the 
Doppler velocities in the region of our interest (cf.. Figs. 1 and 3). 

2.2. Observational results 

Fig. 1 displays the intensity (left), Doppler-velocity (middle), FWHM (right) maps of the 
polar coronal hole where the straight coronal jet has already moved off the limb. The maps 
are related with spectral measurements over 2 pixel x 6 pixel binned data with good signal- 
to-noise ratio. The jet is visible in the left-most part of the intensity map, and it has a 
typical inverted Y-shape structure (cf., within the box in the intensity map). At the 
same place, in the Doppler velocity map, we see the red-shifted base of the jet while the 
plasma moves outward (blue-shift) in the higher plasma column of the jet (cf., within the 
box in the Doppler velocity map). The blue and red-shift measured in the off-limb jet 
corresspond to the line-of-sight component of velocity variations. It seems that the jet is 
slightly tilted towards us. Therefore, at a certain height, the jet plasma which is outflowing 
towards us provides a blue-shift signature, while the downflowing plasma towards the Sun 
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within the jet gives the red-shift. In the FWHM map, we find the increased contribution of 
the line-width in the jet plasma column (cf., within the box in the line-width map). 
It should also be noted that the jet is not visible in Si VII 275.35 line (Log Tf ^ 5.0), 
thereby confirming that it is not formed by cool transition region plasma. The jet is also 
not evident at higher temperatures. It is found that mostly its plasma is maintained at 
the inner coronal temperature of 1.0 MK. Therefore, the jet is best observed in Fe XII 
195.12 A line. We do not conjecture that it is made up of a single temperature plasma. 
In the given observational base-line its plasma temperature can on average be around the 
formation temperature of Fe XII line, i.e.. Log T/ « 6.1. 

We choose a straight slit from the base of the jet in the outward direction along the jet 
plasma column (cf.. Fig. 1). We measure the variation of intensity (top-left panel), Doppler 
velocity (top-right panel), and FWHM (bottom panel) of Fe XII 195.12 A line along this 
chosen path (cf.. Fig. 3). The first position at the bottom of the slit is the base of the jet 
which is referenced as ’zero Mm’. The rest of the slit represents the progressive height along 
the jet scaled in Mm. The total height of the slit along the jet is 56 Mm (cf.. Fig. 3). It 
should be noted that we have considered the estimations of these spectral parameters along 
the jet up to only those heights (w 55 Mm) where spectral fitting is resonable with less 
fitting errors (cf.. Fig. 2). It is clear from these panels that intensity is continously enhanced 
upto 10 Mm and decreases thereafter Doppler velocity changes from its maximum value of 
red-shift about -|-5 kms“^ from the base of the jet to the blue-shift (outflows) at a height 
between 5 — 10 Mm. We find the variation of FWHM along the same path. It is almost 
constant without its much variation upto 10 — 15 Mm (cf.. Fig. 3, bottom panel). Above 
this height, the FWHM increases from 30 mA to 55 A with a significant positive gradient 
(cf., red fitted line in the bottom panel of Fig. 3). Red-line shows third order polynomial 
fit on the FWHM data. 

2.3. Physical interpretations of the observational results 

The small-scale and low-lying bright-point loops form the wide base of the jet. 
They reconnect with pre-existing open field lines and propell jets’ plasma. This 
scenario is in accordance with the model given by Yokoyama & Shibata (1995) about the 
coronal jet formation. The red-shift is maximum (-1-5 kms~^) at 5 Mm height, which means 
the plasma is downfiowing at this location. Above this height between 5 — 10 Mm, this red- 
shift is transformed into blue-shift. This vertical region acts as a reconnection region where 
plasma flows in both directions: the downflows in small-scale loops forming the base of 
the jet and upflows along open field lines of the jet forming its plasma motion higher into 
the corona. This region where downflows (red-shift) convert into outflows (blue-shift) and 
intensity (thus density) tends toward its maximization, indicates that reconnection already 
drove the energy release below of this location near the base of the jet. The reconnection 
propells the plasma outwards along the open field lines of the jet. The maximum observed 
blue-shift is 7.5 kms“^ along the jet. However, it should be noted that the jet is projected 
off the limb at the north-pole. The plasma motions along the jet may have minimum 
contribution to the line-of-sight (LOS) Doppler shift. Therefore, the estimated Doppler 
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velocity should be considered as an apparent lower-bound velocity, which is not as an 
actual speed of the jet. This simply provides the signature of the outflowing plasma along 
the jet. The FWHM also shows the increasing trend beyond 10 — 15 Mm height. It should 
be noted that the FWHM value exactly at the base of the jet at the north polar limb is 
somewhat higher, i.e., 40 mA. We can discard it due to the limb effects. The contribution of 
stray light from continuum near the limb may affect the line profile and its width. Second 
point of estimated FWHM still lies higher around 36 mA at 5 Mm where the highest red- 
shift (downflow) is observed. At 5.0 Mm height, the energy release is likely to be present 
as downflow is maximum at this location. The energy release may cause the line-width 
broadening of the spectral line (Somov 1994). Therefore, the spatial points that lie in the 
lower segment of the jet closer to the energy release site in the spatial range 5—15 Mm 
from Jet’s footpoint, possess comparatively higher values of FWHM. As we move away from 
energy release site at higher heights, the jet’s temperature is uniform near the formation 
temperature of optically thin Fe XH 195.12 A line. Therefore, the contribution to the 
FWHM increment may be mainly due to the increment in non-thermal motions. 

The reconnection generates the outflowing plasma along jet’s open magnetic field lines, 
as well as the downflowing plasma along the small-scale loops forming its base. This region 
(0 — 15 Mm), therefore, may be the most likely place in the lower segment of the jet around 
which the reconnection and related physical processes occurred. Beyond this region, the 
FWHM of the observed spectral line increases significantly due to non-thermal motions as 
evident in the observational base-line. This FWHM broadening may be associated with the 
evolution and growth of the non-linear and impulsively generated Alfven waves propagating 
along the open field lines of the jet away from the reconnection site higher into the corona. 
In the theory as well as observations, Alfven waves are found to be excited in the large-scale 
polar coronal holes (Banerjee et al. 1998, Chmielewski et al. 2013). 

The FWHM also shows varying trend in and around the observed jet. However, these 
trends are not so systematic and closely correlated with intensity and Doppler velocity 
variations with height as we observed within the jet.We do not present these plots. However, 
we have examined the region outside the jet’s plasma column and its inverted Y-shape base 
in its western side. The intensity decreases rapidly as we move outside the limb, which is 
an obvious trend. However, there is no trend in the Doppler velocity. Doppler velocity 
only shows weak downflows. FWHM also does not show any trend in its variation with 
height. In conclusion, present observations demonstrate the signature of impulsive Alfven 
waves along the jet above the energy release site near its footpoint. In the next section, we 
describe our model of impulsive Alfven waves in the 2-D vertical current-sheet to match 
with some observed properties of these wave modes. 

3. Numerical model of impulsive Alfven waves in 2-D vertical Harris 
current-sheet 

Keeping in mind the impulsive generation of Alfven waves due to energy release within 
the observed jet, we model the propagation of a transverse pulse in a 2-D vertical Harris 
current-sheet. Above the reconnection point, perturbations may evolve in the form of Alfven 
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waves propagating at various higher heights as well as the driven coupled plasma flows. We 
consider the height of 5 Mm in our model to launch the reconnection generated velocity 
pulse in a 2-D vertical current-sheet. 


3.1. Governing equations 

In our numerical model of a Harris current-sheet, embedded in a gravitationally stratified 
solar atmosphere, the plasma dynamics is described by the two and half-dimensional (2.5- 
D) time-dependent ideal MHD equations (cf.. Priest, 1982; Chung, 2002): 


Dp V7 

^ = -eV 

(1) 

Du ^ „ 

= -Vp -i-3 X B + gg, 

(2) 


(3) 

De 

- = -leV . V. 

(4) 

< 

to 

II 

0 

(5) 


Here D/Dt = d/dt + v ■ V is the total time derivative, g is the mass density, v is the 
flow velocity, B is the magnetic field, and g = [ 0 , —gQ, 0 ] is the gravitational acceleration, 
that is oriented in the negative y-direction, with gQ = 274 m • s“^. Our model does not 
include the radiative and thermal conductive losses. This model may not fully describe 
the reconnection generated heating and related additional plasma evolution within the jet 
as we do not invoke the non-classical thermodynamical terms (e.g., radiative and thermal 
losses) in the governing energy equation. In the present paper, we do aim to understand 
only the non-linear wave dynamics and associated vertical flows along the jet. 

The current density j in Eq. m is expressed as 

j = —{VxB), ( 6 ) 

go 

where go = 1-26 x 10“® H ■ m“^ is the magnetic permeability of free space. 

The internal energy density, e, in Eq. is given by 


e = 


7-1 

with the adiabatic coefficient which we set and hold fixed as 7 = 5/3. 


(7) 


3.2. Gravitationally stratified Harris current-sheet 

For a static {v = 0) equilibrium, the Lorentz and gravity forces must be balanced by the 
pressure gradient in the entire physical domain 

- Vp-b j X B + gg ^0. ( 8 ) 

The solenoidal condition, V • JB = 0, is identically satisfied by the magnetic flux function, 

A, 


B = V X A. 


(9) 
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For calculating the magnetic field in the vertically oriented Harris current-sheet, we use 
the magnetic flux function A = [0, 0 ,^ 2 ] as (Galsgaard & Roussev, 2002; Jelinek et ah, 
2012) 


A 2 = -BqWcs In 


cosh 1 - 

Wes 


exp - 


(-!)■ 


( 10 ) 


Here the coefficient A denotes the magnetic scale-height. The symbol Bq is used for the 
external magnetic field and Wes is the half-width of the current-sheet. We set and hold fixed 
Wes = 1-0 Mm. Substituting Eq. (10) into Eq. (9) we obtain the equations for the magnetic 
field in the x-y plane as (Galsgaard & Roussev, 2002; Jelmek et ah, 2012) 


Wc 


Bx{x,y) = Ho—In 


By{x,y) = Bq tanh 


cosh 


exp - 


(-!) 


( 11 ) 


exp - 


{-->)- 


( 12 ) 


Equation (8), while written in terms of its components, attains the following form: 
dp{x, y) 


dx 

dp{x, y) 
dy 


+ jzBy{x,y) = 0, 

- jzBx{x, y) + gg{x, y) = 0. 


(13) 


(14) 


Here jz is the only non-zero component of the electric current density j (see Eq. ®), given 
by jz = ^ B)z. The condition of integrability of the above equations leads to 

dg{x, y) 


Mo5- 


dx 


= V ■ (pojzB), 


(15) 


from which we can derive the formulae for the distribution of the mass density (cf., 
Galsgaard & Roussev, 2002; Jelmek et ah, 2012) 
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y 2//0 

gogA [ 
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|.exp(-2|) 


bWc^ 2 
2 goA 2 


'Po- 


cosh 


(17) 


Here g^ and po are integration constants. The corresponding plasma temperature is as¬ 
sumed to be T = 1.25 x 10® K. 


3.3. Initial perturbation 

At the start of the numerical simulation (t = 0 s), the equilibrium is perturbed by the 
Gaussian pulse in the z-component of velocity and has the following form (e.g. Nakariakov 
et al. 2004, 2005): 


Vz{x,y,t = 0) = -Aq ■ - • exp 


+ {y- Lpf 


A 2 


(18) 
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where Aq is the amplitude of the initial pulse, and A = 4 Mm is its width. This pulse 
preferentially triggers Alfven waves. The perturbation point, Lp, is located on the axis of 
the Harris current-sheet, at a distance of 5 Mm from the bottom boundary of the simulation 
region (see full black circle in Fig. 4). We set the detection points, Ld, on the current-sheet 
axis, at the distance between the perturbation and detection points A = |Ld ~ L'p\ = 
10,15, 20, 25, 30, 35,40,45, 50, and 55 Mm. 

Note that in the 2.5-D model, the Alfven waves decouple from magnetoacoustic waves. 
They can be described solely by Vz{x, y, t). As a result, the initial pulse of Eq. (18) triggers 
Alfven waves that are approximately described in the linear case by the wave equation 


d‘^Vz{x,y,t) 2 / ^d‘^yz{x,y,t) 

■ = CA{x,y)- 


(19) 


dt"^ ’ ds^ ’ 

where s is the coordinate along the magnetic line and the Alfven speed, ca, is defined as 


lBl{x,y) + Bl{x,y) 


CA{x,y)=,i , . 

V f^oQ[x,y) 

Here g{x,y) is given by Eq. (16). 


( 20 ) 


3.4. Numerical solutions 

We solve the 2.5-D time-dependent, ideal MHD equations 0-0 numerically, making 
use of the ELASH code (Fryxell et al. 2000; Lee & Deane 2009; Lee 2013). It is now well 
tested, fully modular, parallel, multiphysics, open science, simulation code that implements 
second- and third-order unsplit Godunov solvers with various slope limiters and Riemann 
solvers as well as adaptive mesh refinement (AMR; see e.g. Chung 2002). The Godunov 
solver combines the corner transport upwind method for multi-dimensional integration 
and the constrained transport algorithm for preserving the divergence-free constraint on 
the magnetic field (Lee & Deane 2009). We have used the minmod slope limiter and the 
Riemann solver (e.g. Murawski, 2002; Toro 2006). The main advantage of using AMR 
technique is to refine a numerical grid at steep spatial profiles while keeping a grid coarse 
at the places where fine spatial resolution is not essential. In our case, the AMR strategy 
is based on controlling the numerical errors in a gradient of mass density that leads to 
reduction of the numerical diffusion within the entire simulation region. 

For our numerical simulations, we use a 2-D Eulerian box of its height H = 80 Mm and 
width W = 20 Mm. The spatial resolution of the numerical grid is determined with the 
AMR method. We use the AMR grid with the minimum (maximum) level of the refinement 
blocks set to 3 (7). The whole simulation region is covered by 12884 blocks. Since every 
block consists 8x8 numerical cells, this number of blocks corresponds to 824576 numerical 
cells. 

Note that a spatial grid size has to be less than the typical width of the current-sheet 
along the x-direction and the typical wavelength of the Alfven waves along the y-direction, 
respectively. We find min(Ax) = min(Ay) = 0.04 Mm, which satisfies the above mentioned 
condition (semi-width of the current-sheet is wcs = 1-0 Mm and the estimated minimal 
wavelength is approximately 10.0 Mm). At all boundaries, we fix all plasma quantities 
to their equilibrium values using fixed-in-time boundary conditions, which lead only to 
negligibly small numerical reflections of incident wave signals. 
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4. Results 

4.1. Numerical results 

In a current-sheet, a velocity perturbation is generated most likely due to reconnection 
at a height of 5 Mm (cf., Fig. 4), which evolves with height in the stratified solar atmo¬ 
sphere. The initial pulse excites Alfven waves. These waves, which are associated with the 
perturbations in 14, propagate along the open magnetic field lines to higher altitudes (y- 
direction): see the time evolution of 14 at < = 0,10, 30, and 70 s (Fig. 5). The initial pulse 
(left) spreads into upwardly and downwardly propagating waves which are very well seen 
at t = 10 s (second left). Since these waves follow magnetic field lines, which are diverging 
with height, the left and right profiles of the Alfven waves move apart with y (see the right 
most panel). 

The associated magnetoacoustic waves produce the mass density perturbations at vari¬ 
ous heights and also create the vertical plasma dynamics. Some instantaneous mass density, 
Vy and 14 velocity variations at 15, 35, and 55 Mm heights along the current-sheet above 
the solar surface are shown in Fig. 6. At each height, we have the temporal signatures of 
the perturbations of these parameters due to the evolution of waves. 

4.2. Synthetic line-width increment due to Alfven waves 

In the numerical domain, we have evolution of Alfven waves in the Harris current-sheet, 
which propagate along open magnetic field lines in its stratified atmosphere. These Alfven 
waves are described by the perturbations in the transverse velocity component (14). They 
result in the magneto-plasma as unresolved non-thermal motions (Banerjee et al. 1998; 
Chmielewski et al. 2013). The non-thermal motions due to Alfven waves leave their signa¬ 
ture in the line-profiles at each height in terms of their broadening. The increment in the 
line-width, therefore, provide the signature of Alfven wave amplification (Banerjee et al. 
1998, O’Shea et al. 2005). As mentioned above, the broadening of line-profiles due to Alfven 
waves basically lead to the unresolved motions. This contribution to the spectrum, at any 
particular height in the solar atmosphere, can be due to various wave trains already passed 
through that detection point. However, these motions are not resolved by the spectrometer 
though inherent within the spectrum in terms of line-broadening. Therefore, we average 
over the space (entire pulse width), and time (500 s) of various 14 signals to know the re¬ 
sultant contribution of the unresolved non-thermal motions. The detailed discussions can 
be found in Chmielewski et al. (2013). In this sub-section, we describe the procedure that 
we adopt to convert the simulated velocity signals (14) into synthetic line-width of Fe XB 
195.12 A. At any particular height in the model atmosphere, we average the wave velocity 
amplitude (14) over the entire pulse width between ta — 250 and ta + 250 s, where ta is the 
arrival time of a wave signal to that particular height (Chmielewski et al. 2013). Therefore, 
the averaged transverse motions contribute to the non-thermal unresolved motions at a 
particular height in the model atmosphere as observed in the form of spectral line profile 
variation due to Alfven waves. The averaged wave velocity amplitude is scaled in terms of 
non-thermal speed as = 0.517? by taking into consideration the polarization and the 
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direction of the propagation of Alfven waves with respect to the line-of-sight (Banerjee et 
al. 1998). It should be noted that this scaling can be applicable for longer exposures. Our 
spatio-temporal averaging of velocity signal satisfies this requirement. Using the following 
formula (Mariska, 1992): 


^2 = 


41n2 ( — 


V mi 


+ <jj 


( 21 ) 


we estimate the synthetic line-width (FWHM) of Fe XII (A = 195.12 A; Tj = 1.2 x 10® K) 
as a function of height in the model current-sheet where Alfven waves are present. Here, cr/ 
is the instrumental width associated with the Hinode/EIS spectrometer slit. The observed 
line-width is converted into FWHM and subtracted from instrumental width. While we 
derive the synthetic line-width of Fe XH 195.12 A from numerical data, we ignore this 
factor as it only contains the physical information of thermal and non-thermal motions. It 
is clear from Fig. 7 (diamonds) that the synthetic line width resulting from the contribution 
of Alfven waves grows with height, and its trend matches the variation of the observed line- 
width (cf.. Fig. 2, bottom panel; beyond 10 Mm). 

We conclude that the Alfven waves excited impulsively in our model atmosphere play 
a similar role like the transverse waves excited within the observed polar jet. Therefore, we 
exactly constrain the physical mechanisms, i.e., reconnection generated Alfven waves, which 
also play a role in triggering the observed jet. The detailed comparison is given in the next 
sub-section. The width of the formed line-profile depends both on the thermal and non- 
thermal contributions in the emitted plasma. It is shown that the heating plays a crucial 
role in the formation of spectral lines and its equivalent width (e.g. Peter et al. 2006). In the 
present observational base-line, we have found that our less intense jet does not show much 
temperature variations along its extended part in the corona as observed by Hinode/EIS. 
It is only clearly visible in Fe XH 195.12 A line, while it is not evident in cooler (Si VH 
275.35 A, logT/ = 5.0 MK) as well as hot (Fe XV 284.16 A, logT/ = 6.4 MK) lines. This 
means that the heating episode is quickly subsided after the reconnection generated origin of 
the jet. When it is scanned through Hinode/EIS 2”-slit, it has already reached its maximum 
height in the polar corona with the plasma maintained around inner coronal temperature of 
1.0 MK. Therefore, the thermal width in the formed line profiles of Fe XH 195.12 A (both 
observed and spectral) can be considered as a constant width equivalent to the formation 
temperature of Fe XH line, i.e., Tf = 1.2 x 10® K. The non-thermal contribution and 
resultant line broadening, therefore, result from the Alfven wave amplification along the 
jet in its stratified atmosphere. 


4.3. Comparison with the observational data 

The observed jet is triggered by reconnection between low-lying bipolar loops at its base, 
and pre-existing open field lines of the polar corona. This forms the typical inverted Y- 
shape configuration of the jet, and a vertical Harris current-sheet along the jet body. The 
X-type reconnection point can be created somewhere at coronal heights above the inverted 
Y-shape base of the jet. We model this vertical Harris current-sheet as an upper part 
of the jet where waves and plasma perturbations occur. We assume that the magnetic 
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reconnection generates a velocity pulse which results in Alfven waves propagation along 
the coronal jet. 

In the observations, we find the reconnection point at a height 5 — 10 Mm from the 
base of the jet where maximum downflow speed is converted into blue-shift at higher 
heights along the jet. Up to this height, the measured FWHM is almost constant. It grows 
significantly beyond this height, and maximizes up to 50 Mm height along the jet. As noted 
in Sec. 2, this provides the most likely signature of Alfven wave growth along the jet body 
beyond the reconnection point. 

In our model, the transverse velocity perturbations are generated due to the recon¬ 
nection in the lower part of the Harris current-sheet, which correspond to Alfven waves. 
Although, we include only the reconnection generated velocity pulse in the vertical current- 
sheet of the jet. We obtain the time signatures of Alfven waves at each height, and convert 
it into the unresolved non-thermal velocity. Since mass density decreases with height along 
the open field lines, Alfven waves are amplihed, contributing more to the non-thermal mo¬ 
tions at higher altitudes. This results in the line width growth with height of any observed 
line profile in the stratified solar atmosphere. The synthesized FWHM (line-width) due to 
the contribution of Alfven waves propagating in the current-sheet along the jet body. It 
also shows a similar increasing trend up to y = 55 Mm (cf., Fig. 7). It is constant up to 
y = 15 Mm above the reconnection point. Beyond the reconnection point, the observed 
FWHM profile along the polar jet (Fig. 3, bottom panel) increases. This confirms that the 
Alfven waves are impulsively excited within the typical jet morphology due to magnetic 
reconnection. Alfven waves propagate along the jet’s open magnetic field lines, and further 
leave their signatures in the form of the observed line width variation. Moreover, the ex¬ 
cited slow magnetoaoustic waves perturb the plasma and result in the vertical flows along 
the magnetic field lines, which cause the jet plasma dynamics. In conclusion, our model 
supports the scenario of the evolution of non-linear Alfven waves in the polar coronal jet, 
and associated flows that may contribute to its plasma evolution. However, it should be 
noted that we do not consider the jet’s plasma evolution due to heating or direct j x B 
force that may also contribute initially in its formation. 


5. Discussions and conclusion 

Using the Hinode/EIS data, we observe that magnetic reconnection takes place at the 
coronal height within a polar jet. Above and below the reconnection point, the plasma shows 
upflows (—7 kms~^) and downflows (-1-5 kms~^), respectively. This confirms the formation 
of a typical reconnection-generated coronal jet. We also hnd the line-width increment (Fe 
XH 195.12 A) along the jet body beyond the reconnection point as a likely signature of 
the impulsively generated Alfven waves. The observed line width increment with height 
matches with its theoretically estimated values, obtained by considering the contribution 
of Alfven waves amplitude in the unresolved nonthermal motions of the plasma jet. The 
inference of this matching provides the physics of Alfven waves which are driven within 
the coronal jets. The reconnection within the vertical Harris current-sheet along the jet’s 
body trigger the longitudinal and transverse perturbation, which are associated with fast 
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magnetoacoustic and Alfven waves. The Alfven waves leave their signatures collectively 
in the form of spectral line broadening with height, while the fast magnetoacoustic waves 
form the plasma flows in the jet. Linear and non-linear Alfven waves can also be generated 
at the photosphere by granular motions. The region between photosphere and transition 
region acts as an Alfvenic resonator (e.g., Matsumoto & Shibata 2010). In the present 
case, we do not find any observational signature of Alfven waves below 10 Mm within 
the jet. Therefore, we conjecture that the Alfven waves are impulsively triggered above the 
reconnection site lying between 5 — 10 Mm above the footpoint of the jet. The photospheric 
Alfven waves generated by granular motions have typical initial amplitude of a few kms~^. 
When they propagate higher into the stratified atmosphere, they do amplify. However, 
in the present case, the Alfven waves are triggered in situ at coronal height and possess 
the larger instantaneous amplitudes (> 10 kms“^). Therefore, we abbreviate them as 
impulsively generated non-linear Alfven waves propagating along the polar jet. 

The underlying reconnection can drive Alfven waves which excite the vertical plasma 
flow. We did not aim to understand the temporal evolution of jet plasma and its heating as 
previously explored by various models (e.g. Moreno-Insertis et al. 2013). Instead, we aim to 
understand the signature of Alfven wave dynamics and associated vertical flows along the 
jet. We find the clues of both respectively in the form of line-width increment and Doppler 
velocity distribution (blue-shift) along the jet at various heights. For this purpose, we 
select the Hinode/EIS scanning spectral observations to understand the spatial variation of 
emission (flux) and plasma properties (flows and width) along the jet. The jet was initially 
launched due to reconnection which may also liberate the heat. The heat flux may be 
subsided during the whole half-life time of the jet when it has already reached a maximum 
altitude in the polar corona. The impulsive launch of the jet might leave other physical 
processes for rest of the time in space along with the jet, i.e., Alfven waves and associated 
vertical flows. We have observed these two physical processes through EIS spectroscopic 
scan data, which also confirm the model. Our model is simple in the sense that it does 
not include the radiative or thermal conductive losses. As a result, this model may not 
fully describe the reconnection generated heating and related additional plasma evolution 
within the jet as we do not include the thermodynamical losses in the governing energy 
equation. In general, the non-consideration of these terms may not affect the properties 
of Alfven waves that evolved along the jet. However, they do affect the mass density and 
temperature evolution of the jet plasma. In our case, the jet is visible in the inner coronal 
temperature of 1.0 MK, and it is not a very hot evolution of the jet. Moreover, weak EUV 
emissions also suggest that plasma is not well evolved along the jet spine and it is not a 
bulky jet likewise X-ray jets and surges. In order to understand and compare the properties 
of the triggered Alfven waves, we consider the simplihed ideal situation in our model by 
ignoring the losses. These thermodynamical losses will be considered in our future works 
in which we will examine their effects on the evolution of temperature and density within 
the jet. 

There exists only one report in the imaging X-ray observations from Hinode/XRT, 
which shows that the small amplitude Alfven waves propagating from photosphere into 
the corona can drive the X-ray jets (Cirtain et al. 2007). Such important detections are 
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debated by Van Doorsselaere et al. (2008a, 2008b). There are several attempts to model 
the polar jets in terms of the evolution of Alfven waves (Kaghashvili 2009; Pariat et al. 
2009, 2015; Hollweg & Kaghashvili 2012, and references cited therein). These models are 
purely analytical/numerical without any observational support. Nevertheless, they typically 
reveal the physics of the jet forming regions as to how such waves can be excited within the 
typical magnetic field configurations of these jets. Nishizuka et al. (2008) have proposed 
a model of the coronal jets by extending the model of Yokoyama & Shibata (1995), and 
found that direct reconnection generated forces can drive both the plasma jet and Alfven 
waves at its base. Recently, Kayshap et al. (2013) reported that the reconnection in the 
lower atmosphere can generate the pulse train and associated slow shocks to drive the 
surge plasma. Shibata et al. (2007) have also suggested three types of the processes for 
the jet formation caused by reconnection at different heights. Depending upon the height 
of magnetic reconnection in the pre-existing magnetic field configuration of the jet, the 
hot X-ray, EUV, and cool jets as well as Hq, surges can be triggered respectively in the 
coronal, transition region, and chromospheric reconnection. Therefore, reconnection can 
be an efficient mechanism to drive the waves/pulses within the confined magneto-plasma 
system (i.e., various jets) which can also play a role to power them. Alfven waves are well 
studied both in theory and observations in driving the small-scale chromospheric jets (e.g., 
spicules; Kudoh & Shibata 1999; McIntosh et al. 2011). 

In the present work, we have specifically made an effort to find the signatures of recon- 
nction within a jet at a certain height from where the impulsively excited Alfven waves 
and plasma motions are evident. We find the matching of the observed scenario with the 
given physical model of impulsive excitation of Alfven waves within the vertical Harris 
current-sheet. Our model has the appropriate temperature conditions as well as stratified 
atmosphere in 2.5-D framework along with a consideration of the transition region which 
plays an important role in reflecting the Alfven waves. The transition region has an impor¬ 
tant implication on contribution in non-thermal motions of the corona and in the formation 
of spectral lines provided they are excited below the transition region. In the present case, 
the site for impulsive energy release lies in the inner corona. Therefore, the excited Alfven 
waves propagating in the outward direction will not have any influence on the transition 
region. Moreover, the appropriate temperature conditions also set the stratification within 
the model atmosphere, which will result in the appropriate growth of the Alfven waves. 
Therefore, its contribution to the non-thermal motions and corressponding synthetic line 
broadening are adequately examined and studied. They also match with the observed line 
width variation along the polar jet. 

In conclusion, the numerical simulations show that the perturbations in the transversal 
component of velocity above the reconnection point within the current-sheet (the upper 
part of the jet and its reconnection site) can excite Alfven waves impulsively. This can also 
trigger the plasma jet higher in the polar coronal hole as vertical plasma flows associated 
with magnetoacoustic waves. The inferred contribution of simulated Alfven waves and their 
synthetic line widths match with the observed non-thermal line broadening along the jet. 
This provides a direct spectroscopic signature of the reconnection generated Alfven waves 
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within the polar jet in its typical magnetic field configuration (i.e., the Harris current-sheet 
and its X-type reconnection site). 
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Fig. 2. Gaussian fit to the spectra derived on the chosen path along the jet. The white slit 
on the jet (cf., Fig. 1) is divided in 14 super pixels each of the size of 2 x 6 pixel square. 
Line profile is derived from each super pixel (dotted-blue), and corressponding estimation 
of plasma parameters are made by fitting the Double Gaussian (solid-blue). 
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Fig. 3. The variation of intensity, Doppler velocity, and FWHM of Fe XII 195.12 A line 
along the path chosen on the polar jet from its base outward. Red-line shows third order 
polynomial fit to the FWHM data. 
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Fig. 4. The mass density distribution in the Harris current-sheet, where the black contours 
represent the structure of magnetic field. On the right-hand side of the figure both (a) 
vertical (for x = 0 Mm) and (b) horizontal (for y = 37.65 Mm, which is magnetic scale 
height) 

slices in mass density are shown, as well as the details of the part of the simulation region illustrat¬ 
ing the computational grid, with Adaptive Mesh Refinement (AMR). At the altitude y = 5 Mm 
the perturbation point, Lp, is marked by full black circle. 
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Fig. 5. The time evolution of velocity component 14, corresponding to Alfven waves, at 
t = 0, 10, 30, and 70 s from left to right. 
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Fig. 6. The top, middle and bottom panels display respectively the spatio-temporal aver¬ 
aged mass density g, velocity components Vy and 14 at three different heights 15,35 and 
55 Mm above the solar surface. 



Fig. 7. The variation of synthetic line width with height (diamonds) of Fe XII 195.12 A as 
derived by considering the contribution of propagating Alfven waves in the Harris current- 
sheet, and the observed line-width broadening along the polar jet (asterisks). 
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